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ABSTRACT. We calculated models of massive AGB stars with a self-consistent coupling of time- 
dependent mixing and nuclear burning for 30 isotopes and 74 reactions. Overshoot with an exponentially 
declining velocity field was considered and applied during all stages of evolution and in all convective 
regions. Very efficient 3"^ dredge-up was found even overcompensating the growth of the hydrogen- 
exhausted core after a few thermal pulses. Hot bottom burning occurs for M > 4Mq within the 
sequences with overshoot. Carbon star formation in these more massive AGB stars is delayed or even 
prevented by hot bottom burning despite the very efficient dredge-up. With the simultaneous treatment 
of mixing and burning the formation of Li-rich AGB stars due to the Cameron-Fowler mechanism was 
followed. For a 6Mq model the maximum Li abundance was found to be e(^Li) ^ 4.4. 

1. Introduction 

In recent stellar evolution calculations Herwig et al. (1997) have considered diffusive 
overshoot for all convective boundaries which leads to a considerable change in the 
models. This method provides for Asymptotic Giant Branch (AGB) stars a sufficient 
amount of dredge-up to form low-mass carbon stars as required from observations. 
Additionally it leads to the formation of ^'^C as a neutron source to drive the s-process 
in these stars. More details are given in Herwig et al. (1999). 

However the consequences for massive AGB stars, which suffer from hot bottom 
burning, remain to be investigated in detail. Hot bottom burning, i.e. the penetration 
of the convective envelope into the hydrogen burning shell, provides lithium-rich AGB 
stars and may delay or prevent the carbon star stage by turning ^^C into ^^C and 
-'^■*N. We have extended our computations of AGB stars with diffusive overshooting to 
hot bottom burning models including a self-consistent coupling of nucleosynthesis and 
mixing. 



2. Computational details 

We used the evolutionary code of Blocker (1995). Nuclear burning was accounted for by 
a nuclear network covering 30 isotopes and 74 reactions up to carbon burning. We con- 
sidered the opacities of Iglesias & Rogers (1996) supplemented with those of Alexander 
& Ferguson (1994) for the low temperature regime. 



Convection was treated within the mixing length theory (Bohn-Vitcnsc 1958) with 
a mixing- length parameter a = 1.7. Mixing of chemical elements was treated by solving 
a diffusion equation. Overshoot was taken into account according to the prescription of 
Herwig ct al. (1997) which is based on the hydrodynamical calculations of Frcytag ct al. 
(1996). Freytag et al. (1996) showed that mixing takes place well beyond the classical 
Schwarzschild border due to overshooting convective elements with an exponentially 
declining velocity field. In the overshoot region the corresponding diffusion coefficient is 
given by Dqs = VQ-Hp-exp with vq: velocity of the convective elements immediately 
before the Schwarzschild border; z: distance from the edge of the convective zone; /: 
the overshoot efficiency parameter. We used an efficiency parameter of / = 0.016 as 
appropriate to match the observed width of the main sequence. Diffusive overshooting 
was considered in all convective regions during the complete evolution. 

Abundance changes have been treated self-consistently by coupling time-dependent 
mixing and nuclear burning of all chemical elements, i.e. by solving 



The simultaneous treatment of burning and mixing is essential to follow, e.g. the ''Li- 
production in luminous AGB stars via the Cameron-Fowler mechanism (Cameron & 
Fowler 1971, sec also Sackmann & Boothroyd 1992). The sclf-consistcnt solution of time- 
dependent burning and mixing processes is, however, associated with a considerable in- 
crease of computing time due to the size of the problem (30 isotopes, ~ 2000 mass shells). 
We calculated sequences for initial masses of SM© to 6M0 and {X,Y) = (0.70,0.28) 
from the pre-main sequence stage up to the AGB and through the thermal pulses. In 
order to disentangle the influence of mass loss from the one of overshoot we refer in this 
study only to sequences with small mass-loss rates (Rcimers 1975, f] = 1.0). Complete 
sequences will also consider stronger (and better suited) rates (cf. Blocker 1995). 

3. Thermal pulses and third dredge up 

On the upper AGB the helium burning shell becomes recurrently unstable raising the 
thermal pulses (Schwarzschild & Harm 1965, Weigert 1966). During these instabilities 
the luminosity of the He shell increases rapidly for a short time of 100 yr to 10^ to 10^ Lq. 
The huge amount of energy produced forces the development of a pulse-driven convection 
zone which mixes products of He burning, carbon and oxygen, into the intershell region. 
Because the H shell is pushed concomitantly into cooler domains H burning ceases 
temporarily allowing the envelope convection to proceed downwards after the pulse, to 
penetrate those intershell regions formerly enriched with carbon and to mix this material 
to the surface (3'''^ dredge up, see, e.g.. Blocker 1999 for a recent review). Overshoot leads 
to an enlargement of the pulse-driven convection zone and to enhanced mixing of core 
material from deep layers below the He shell to the intershell zone ("intershell dredge- 
up" ) and to a deepening of the envelope convection. If the determination of convective 
boundaries is solely based on the Schwarzschild criterion as in our case, dredge up 
can easily be obtained if some envelope overshoot is present to overcome the H/He 
discontinuity. The total amount of dredge up, however, depends mainly on the strength 
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Fig. 1. Temporal evolution of the extension of convective regions for a 6M0 model during the 
second thermal pulse. The hatched regions refer to the bottom of envelope convection and to 
the pulse-driven convective zone of the He shell. The dashed-dotted line refers to the core mass 
Mh, and the dashed one to regions of maximum energy production enc.max in the He shell. 



of the former intershell dredge-up (cf. Herwig et al. 1999). As in the case of low-mass stars 
intershell dredge-up leads to considerable changes of the intershell abundances. After 
intershell dredge-up the abundances (mass fractions) of He, C, O amount to (40,40,14) 
instead to (70,25,2) as in non-overshoot sequences. 

Fig. shows how the lower boundary of the envelope convection as well as the pulse- 
driven convection zone evolve during a pulse for a 6 Mq model. The introduction of 
overshoot and its application to all convective boundaries provides 3'^'^ dredge up even 
for low-mass AGB stars (Herwig et al. 1997, 1999) and we found very efficient dredge up 
for higher masses as well. With the dredge-up parameter A being the ratio of material 
mixed up during the dredge up and the material burnt during two consecutive thermal 
pulses we find A > 1 already during the very first pulses. Correspondingly, one observes 
a decrease of the mass of the hydrogen-exhausted core Mh (see Fig. ^) . 

Due to the overshoot we found high temperatures at the bottom of the pulse-driven 
convection zone. As demonstrated in Fig. |2| one obtains for a 3M0 sequence temper- 
atures already close to 300 • 10^ K, the threshold temperature of the neutron source 
^^Ne(a,n)^^Mg to operate. More massive models show even temperatures of 350 • 10^ K 
and more after a few pulses. Thus, the s-process nucleosynthesis in these stars will be 
governed by both the ^'^C(a,n)^^0 and the ^^Ne(a,n)^^Mg neutron source. 
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Fig. 2. Temperature vs. mass for convective regions for different initial masses. The hatched re- 
gions refer to the bottom of envelope convection and the "spikes" to the pulse-driven convective 
zone of the He shell. 



4. Hot Bottom Burning 

In more massive AGB stars (il/initiai ^5 4M0) the convective envelope becomes so ex- 
tended downwards that it can cut into the hydrogen burning shell during the interpulse 
phase (hot bottom burning (HBB) or envelope burning; Iben 1975, Scalo et al. 1975). 
Temperatures in excess of 50- 106 K (see Fig. |) are reached at the base of the convective 
envelope and material burnt there is immediately mixed to the surface. HBB models do 
not obey Paczynski's (1970) classical core-mass luminosity relation but, instead, evolve 
rapidly to very high luminosities (Blocker & Schonberner 1991). Due to CNO cycling 
of the envelope ^^C can be transformed into ^'^C and ^^N. Consequently, a low ^^C/^'^C 
ratio is a typical signature for HBB which can prevent AGB stars from becoming carbon 
stars (Iben 1975; Renzini & VoU 1981; Boothroyd et al. 1993, Frost et al. 1998). 

The evolution of different CNO isotopic ratios at the surface of our 5 Mq sequence 
is shown in the lower panel of Fig. bI The ^^C/^'^C ratio reaches its equilibrium value 
(w 3) after the first pulses and ^^'Np-'^C starts to increase whereas ^^C/^^0 decreases. 
Although dredge-up is very effective leading even to a core-mass decrease, HBB is so 
efficient that it prevents or delays the star to become a carbon star. 

The evolution of the surface abundances depends on the competition between HBB 
and dredge up. Both effects depend on the envelope mass, i.e. on mass loss. Consequently, 
the C/0 ratio reached at the end of the AGB evolution depends crucially on the question 
which process will shut down first. The critical minimum enevelope mass for efficient 
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Fig. 3. Temporal evolution of the core mass Mh, the shell luminosities Lh and Lhb, and surface 
CNO isotopic ratios (from top to bottom) for a 5 M0 sequence. 

HBB appears to be larger than the one necessary for dredge up (Frost et al. 1998), and it 
depends on the remaining dredge-up eiEciency if finally carbon stars can be formed. Our 
calculations show that dredge up even operates during the post-AGB stage if overshoot 
is considered with important conclusions for the formation of hydrogen-deficient stars. 

Another direct consequence of HBB is the formation of lithium-rich stars (Scalo et 
al. 1975) via the Cameron & Fowler (1971) mechanism: The timescale of the /3-decay of 
^Be which is produced at the bottom of the envelope convection is comparable to the 
convective timescale. Therefore, before ''Be is burnt in the lower part of the envelope 
it can be mixed up to cooler layers where it decays to produce ^Li which, in turn, is 
then convected to the surface. This mechanism is effective for temperatures between 30 
and 80 • 10^ K at the bottom of the envelope. The production of Li-rich AGB stars can 
only be followed with a simultaneous treatment of mixing and burning. Otherwise, ''Be 
would be burnt at the bottom of the envelope before it can be mixed up. 

Since the luminosity is a unique function of the base temperature and Li is only pro- 
duced in a certain temperature range, Li-rich stars are only found between Mboi ~ -6 to -7 



(Sackmann & Boothroyd 1992) in good agreement with observations of the Magellanic 
Clouds (Smith & Lambert 1990, Plez et al. 1993). 

On the other hand, the majority of Li-rich galactic AGB stars are observed at much 
lower luminosities, Mboi ~ —3.5 to —6 (Abia & Isern 1997). They are, however, carbon 
stars. Since ^Li enrichment requires Tbottom ^ 30 • 10^ K whereas ^^C is efficiently 
destroyed via CNO cycling for Tbottom k, 70 • 10^ K one finds in standard evolution 
calculations only a narrow mass, age and brightness range where AGB stars should 
be both lithium and carbon rich indicating that possibly additional mixing processes 
operate in these stars (Abia & Isern 1997). 

Recently, Ventura et al. (1999) calculated AGB models for the LMC with M < 4M0 
considering overshoot in a similar manner as Herwig et al. (1997). They found Li-rich 
carbon stars only for M = 3.5 and 3.8 Mq, i.e. for the upper part of the observed 
luminosity range. Our 4Mq model {Z = 0.02) becomes a carbon star after 10 thermal 
pulses due to very efficient dredge up and can be expected to become lithium rich 
during the further course of evolution due to the sufficient increase of envelope base- 
temperatures. The 3 M0 model becomes a carbon star after a couple of thermal pulses 
as well but seems to miss the Li-rich stage. These findings are in line with those of 
Ventura et al. (1999). Thus, overshoot alone appears not to be able to explain galactic 
Li-rich carbon for the whole brightness range observed. 

Fig. ^ illustrates the considerable ^Li enrichment during the first thermal pulses for 
the 6 Mq sequence. Within the first 3 pulses the Li abundance increases by more than 5 
orders of magnitude to reach almost 10 times the main sequence value. The maximum 
lithium abundance is e(^Li) — log [n('^Li)/n(H)] +12 « 4.4 in agreement with the results 
of Sackmann & Boothroyd (1992). 

Mazzitelli et al. (1999) calculated lithium production in AGB stars with application 
of overshoot to all convective regions as well. The overshoot treatment is comparable to 
the method of Herwig et al. (1997) but convection is dealt with the Canuto & Mazzitelli 
(1992) prescription. They obtained a similar evolution of lithium enrichment and base 
temperatures but found for a solar metallicity 6 Mq model a complete penetration of 
the hydrogen burning shell by the convective envelope resulting in a lack of helium 
accumulation in the intershell region and corresponding prevention of thermal pulses. A 
corresponding quenching of thermal pulses was not found for the present 6 Mq sequence 
(see Fig. I). 

5. Conclusions 

We calculated models of massive AGB stars (solar metallicity) with a self-consistent 
coupling of time-dependent mixing and nuclear burning for 30 isotopes and 74 reac- 
tions. Treating exponential overshoot (Freytag et al. 1996) as in Herwig et al. (1997) 
and applying it to all convective regions we find very efficient 3"^ dredge-up with A > 1 
after a few thermal pulses. Overshoot leads to intershell dredge- up and provides con- 
siderably changes of the intershell abundances strongly enriched in carbon and oxygen, 
i.e. [He,C,O] = [40,40,14] by mass. These intershell abundances determine the strength 
of the third dredge up, whereas envelope overshoot is mainly required to penetrate the 
H/He discontinuity. Temperatures at the bottom of the flash-driven convection zone of 
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Fig. 4. Temperature at the base of the convective envelope for models with initial masses of 4, 
5 and 6Mq The strong increase of Tbottom indicates the efficiency of HBB. t = refers to the 
I/He maximum of the first thermal pulse. 
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Fig. 5. Upper panel: Evolution of the shell luminosities Lh and Luc, and of the surface 
abundance of ^Li with respect to its initial value (^LiinUiai = 9.35-10"^) during the first thermal 
pulses of the M = 6Mq sequence. Lower panel: Corresponding evolution of the temperature at 
the bottom of the envelope convection. 



the helium burning shell are in excess of 300 • 10® K after a few pulses activating the 
^^Ne(a,n)^^Mg reaction for s-process nucleosynhesis. 

Hot bottom burning occurs for M > 4M0 within our sequences with overshoot. 
Only the 4Mq model becomes a carbon star whereas carbon star formation in more 
massive AGB stars is delayed or even prevented by hot bottom burning despite of the 
very efficient dredge- up. With the simultaneous treatment of mixing and burning the 
formation of Li-rich AGB stars due to the Cameron-Fowler mechanism was studied. 
Already the 4 model becomes lithium-rich leading to the formation of a lithium-rich 
carbon star. For M > 6 M© the maximum Li abundance was found to be e(^Li) ps 4.4. 
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